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VIII. — On the Equilibrium and Motion of Solid and Fluid Bodies. By the Rev. 

Samuel Haughton, F. T. C. D. 



Bead May 25, 1846. 



A HE object of the present paper is to deduce from simple physical considerations 
the laws of equilibrium and motion of elastic solid and fluid bodies, by the 
method followed in the Mecanique Analytique of Lagrange, which possesses the 
remarkable advantage of giving, by the same analysis, the general equations of 
any system, together with the particular conditions to be fulfilled at the limits. 
This method is particularly valuable in such cases as the present, where the 
problem is, to determine the conditions of equilibrium or motion of an indefinite 
number of material points, situated indefinitely near each other, and acting 
according to certain laws ; although it may be admitted that the ordinary methods 
usual in mechanics are preferable to Lagrange's method, in cases where the 
points of application of the forces of the system are definite in number, and 
situated at finite distances from each other. The first application of Lagrange's 
method to the problem of elastic solids was made by Navier (Mem. de FInstitut. 
Tom. VII.), who discussed the laws of equilibrium of a homogeneous uncrystal- 
line solid, but does not appear to have ever undertaken the general problem, or 
to have obtained any dynamical results. The present paper is an attempt to 
apply the same method of Lagrange to the general case of material substances, 
whether fluid or solid, homogeneous or heterogeneous, and whether possessed of 
a crystalline structure or not ; and more particularly to investigate the general 
dynamical laws of solid elastic bodies, and the conditions at the limiting surfaces 
which bound the solid. 

In the method of Lagrange the forces in action in such a system as a con- 
tinuous body must be divided into two parts — the external and the internal forces 
of the system ; and the general equation of equilibrium is determined by express- 
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152 The Rev. S. Haughton on the Equilibrium and 

ing that the sum of the moments (as defined by Lagrange) of the external forces 
is equal to the sum of the moments of the internal forces ; hence the general 
equation of the system is 

SSS (xB£ + tfn + zS£)dm = SSS Zvdxdydz, (1) 

where x,y, z are the coordinates of the position of rest of any molecule (to), and 
x -f- £, y + 1), z -f- £ , are the coordinates of the same molecule when displaced ; 
and the function v depends on the medium in such a manner that 8v is the sum 
of the moments of internal forces acting at any point of the medium ; and 
SSS fadxdydz represents the sum of all the moments arising from internal forces. 
It appears, therefore, that in investigating the laws of any medium, it is necessary 
to know the form of the function v ; being given the form of v, the laws of the 
medium and of the actions of its molecules on each other are known, and vice 
versa, if we know the nature of the medium, we can determine the form of v. 
Either of two methods may, therefore, be followed, in discussing the conditions 
of a system such as a solid or fluid body ; the form of the function v may be 
assumed, and the truth of the assumption afterwards verified, — or the form of the 
function v may be determined from whatever assumption we make, as to the mode 
of action of the molecules of the system on each other : the latter method is used in 
this paper. The assumed definition of solid and fluid bodies is the following : that 
they consist of molecules, which act on each other only in the line joining them; 
i. e. the whole of the action of one molecule on another is a force in the 
line joining the molecules, and none at all in a direction perpendicular to 
that line. This force is also assumed to be a function of the distance of 
the molecules, and in the case of crystalline structure also of the direction of 
the line joining them. All solid and fluid bodies are either homogeneous or 
heterogeneous in their structure : in homogeneous bodies the molecular forces 
will, ceteris paribus, be the same at all points of the body, but in heterogeneous 
bodies they will vary from point to point of the body, and become functions 
of (x, y, z). A crystalline structure is different, and consists in this, that 
the action of the molecules surrounding any molecule (to), depends upon the 
direction in which they lie, and not merely on their distance from (to). In order 
to avoid misconception, it is necessary to state that the molecular forces alone are 
included in the internal forces ; the mutual gravitations of the particles towards 
each other, although they might with propriety be called internal forces, are 
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here included among the external forces, as they are quite different from the 
molecular forces, and act at finite distances ; and any statement made in reference 
to the function "V is to be understood only of the molecular forces, which alone 
are included in it. Also, all such forces as those arising from the pressures of 
the sides of vessels containing a fluid, &c. are excluded from the function v. 

Let (p, a, /9, 7) be the original length and direction of the line (mm') joining 
two molecules ; suppose these to become (/>-f- />', a-f- <*', /3 -f- ft, 7 + 7'), where 
(/>', a',ff,<y') are the changes in length and direction ; it is evident that the force 
acting between m and m' will be, in general, a function of (p, a, /3, 7, p', a', p y y), 
and also of (ar, y, z) ; or, if (/) denote this force, it may be represented by 

/=F + 2F/ + 3F 2 / 2 -f-&C. 
' A,*' 4- A 2 a' 2 -f &C. 

+ J b/ + b 2 /3" + &c I (2) 

. «*/ + c 2 7' 2 + &<-•. J 
f , f,, &c, A„ b„ c,, &c, being functions of (p, a, /3, 7, x, y, z) in general ; and 
the part of (f) depending on (/»') arises from the change of distance ; while the 
part depending on (a', /S', y) arises from the change of direction. If we con- 
sider, however, that the change of force produced by a change in direction is 
continuoics, while an indefinitely small change of distance developes a finite 
molecular force, it will be evident that we may neglect the part of (/) depending 
on change of direction, together with the part arising from high powers of (//) ; 

consequently 

/=f + 2f i/ / + &c, (3) 

where f„, f„, &c, are functions of (/», a, /3, 7), and implicitly of (x, y, z), if exter- 
nal forces act on the system, or if the body be naturally heterogeneous. 

As the force (/) is a force tending to alter p-\-p', its moment will be f.tp ; 
hence, neglecting the higher powers of //, we shall have the sum of moments of 
internal forces at a given point, 

av = S(F 8/-f.F 1 8( P ') 2 ); 
and, finally, 

v = TSS'M + F ^ ?™<>dpd°d*> (4) 

and 

v = v + v, ; 

v arising from the terms f , and v l arising from the terms f,. 

y 2 
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This value of the function v contains all solid and fluid bodies, together with 
the class of bodies which cannot be said to be either solid or perfectly fluid, and 
which are generally called «' viscous" fluids. If we knew sufficiently well the 
difference of molecular structure of these different classes of bodies, we should 
be able to determine the proper v for each of them ; but in the absence of accu- 
rate knowledge on this point, we must only give the best definition we can of 
solids and fluids, and the agreement of the secondary laws deduced from the defi- 
nition, with the observed laws of these bodies, will be sufficient proof of the cor- 
rectness of the definition founded on their molecular structure. 

I define a solid body to be of such a nature that the part of the molecular force 
depending on F does not exist, so that in a solid body, f and, consequently, v , 
is zero. In all other bodies, liquid, gaseous, and viscous, the function v will 
contain a part v w which does not enter into the function peculiar to solid bodies. 
The definition of a solid by means of the assumption that f and v are zero, is the 
same thing as supposing that in a solid body the molecular forces equilibrate each 
other, without the aid of external forces ; while in fluid bodies the molecular 
forces do not equilibrate themselves without the aid of external forces and pres- 
sures ; and this notion of a solid seems to be correct from the nature of the mole- 
cular forces themselves, which appear to be attractive up to a certain distance ; 
but if the distance be still further diminished, the force becomes repulsive. The 
condition v„ = 0, may then arise from the molecules of a solid body being placed 
at such distances from each other that they are held in equilibrium by the action 
of powerful, attractive, and repulsive forces acting in opposite directions. How- 
ever, our ignorance of the internal structure of bodies is so great, that such 
considerations as these would be, by themselves, an unsafe foundation for a theory 
of solids and fluids ; and the best course is, to use them only as indications of 
the correct definitions, which must themselves be ascertained by the agreement 
of the results they lead to, with the known laws of solid and fluid bodies. 

In following this method in the present paper, I hope to show that the 
presence of v will produce the equations of hydrostatics and hydrodynamics, 
while the term v, will give all the laws of solid bodies ; so that we have the fol- 
lowing functions; — for solids, v = v,; 

for perfect fluids, v = v ; 

for viscous fluids, v = v -f- v,. 
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In perfect fluids the terms arising from v, may be neglected in comparison 
with those produced by v„» at least in hydrostatics ; though I think it probable 
that in hydrodynamics the use of the function v -}- V, would be sometimes neces- 
sary. Also in viscous fluids v , and the terms produced by it, will be very small, 
so that in this case we should always be compelled to use the function v -{- v,. 

I now proceed to the investigation of the equations of equilibrium and motion 
of fluids, or bodies in which v = v . The general equation of equilibrium and 
motion in this case will be 

SSS ( xS * + *^ + *8£) An = SSS Sv dxdydz, (5) 

where 

The value of p, in terms of £, % £ , &c, to be substituted in this formula, is thus 
found. Let x,y,z, x + a, y -j- b, z -J- c, be the coordinates of the positions of 
rest of m and to' ; then in the altered position of the molecules, if x,y, z become 
x -\- £, y -\- r\, z + £ , the coordinates of to' will become 

.+«+.+£. + !»+§• 

and consequently p -\- p', being equal to the square root of the sum of squares of 
differences of coordinates of its extreme points, will be expressed by the formula, 

f W{(-+£+3»+*)+(»+i-+$»+i-)+(-**3»IW 

from which, by neglecting the smaller quantities, and assuming 

_ dri d£ __ d£ d% _ d% dq 
V ~Tz + dy' y ~dx + dz' w -dy + dx' 
we obtain 

p-\-p' =p + p(-^- C0S 8 O-f--^C0S S j8-f--^ cos 2 7+ ucos/3cos7-f-vcosacos7-f-wcosacos/3 ), 
and, finally, 
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p' = pi T-COS 2 a -f" -T- COS 2 /3 -\- -^•COS 2 7+ UCOS/3COS7 + VCOSOCOS74- wcosacos/3 j. (6) 

Substituting this value of p' in (5), we obtain 

ucos/3cosy -j- v cosa cos? -J- wcosacos/8 j p 3 sm0dpd8d<p ; 

or, since, 

cos* = cos^ sin0, cos/3 = sin^ sinfl, COS7 = co&O, 

i. e. v is a homogeneous function of the first order of the six quantities, 

dt- dv\ dg 

dx dy dz 
where 

Aj = SSS f . cos 2 0sin 3 p 3 dpd6d<p> D, = SSS f . sirupsitfecose p 3 dpd6d<p, 

b, = SSS F o .sin 2 0sin 3 p 3 dpd6d<p, e, = SSS r o .cos<psin l 0co<i6 p 3 dpd8d<p, 

c, = SSS v .cos' i 6sm6.p 3 dpd6d<p, f, = SSS F o .sin0cos0.sin 3 tf/> 3 d/»d0cty>; 

but, since in fluids, the structure is not crystalline, F will not be a function of 
(0, tp), but only of p; and by integrating twice with respect to 6 and <p, we 

obtain, 

Dj = e, = f, == 0, 

A, = B t = c, —p = -£ ^ F p 3 d/> ; 
hence we deduce, finally, for the value of v , the expression 

^(§ + g+f> (7 > 

Substituting this expression for v ffl in the general equation of fluids (5), we obtain 
by the calculus of variations, the following result, 

SSS (xB? + yhrj + zS£)dm = §pl&ydz + §pfydxdz + \\pl£dxdy 

-fiS^+^+i*)*** < 8 > 
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the triple integrals giving the equations of equilibrium, and the double integrals 
giving the conditions at the limits. If the density be expressed by e, so that 
dm = edxdydz, we have 

_ €X -^ _ £Y -^ -ez~& «1) 

€X ~dx' dy> eZ ~dz' W> 

which are the well-known equations of hydrostatics. It appears from these equa- 
tions that if the forces (x, y, z) be zero at all points of the fluid, the quantity p 
must be constant ; and vice versa, ifp be constant for all points of the fluid, the 
forces x, y, z must be zero ; in such a case the function v will not give any inde- 
finite equations of equilibrium, but only the particular condition at the limits ex- 
pressed by the double integrals, which is the same as that found by Lagrange, 
and expresses that there must be a constant normal pressure at each point of the 
limiting surface, equal to p, in order that the fluid should remain in equilibrium. 
The hydrodynamical equations corresponding to (9) are 

e^-ex-4-^ 6 ^- eY -+.^ e £l- ez -L.± no) 

which are the equations commonly used in hydrodynamics, and, together with the 
equations (9), are usually deduced from the principle of equality of pressure, of 
which principle they are merely the mathematical expression : but the fact that 
this principle is often not true in the case of fluids in motion, would lead us to 
suppose that though the equations may do well enough as a first approximation, 
yet that in hydrodynamics they are insufficient ; in fact we have, in general, 

v = v + v, -f v 2 + &c. 

and the terms v,,, v„ &c. give rise, in the equations of equilibrium and motion, to 
differential coefficients of the first, second, &c, order ; and in case the supposition 
v = v becomes insufficient, we must then make v = v + v„ and add to the 
differential equations of equilibrium and motion, arising from v„, the additional 
terms produced by v,. I shall return to this subject in a subsequent part of this 
paper, when I have determined the form of the terms produced by v,. 
It appears from (4) that 

v,=s;s;s> I (/' / )v s sin^^#, 
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from which the following value for v, may be deduced, by substituting for (/>') its 
value taken from (6), 

»• = i (l) , +<|)'+<D'+- , + "*■+ «■ <») 
+ 2 !<4+4+-l)+<4+4 + 4)+"(4+4 + 4)S 

and if we make dw zz p 4 sm0dpd9d<p, we shall have the following values for the 
coefficients of the function v x : 

a = 2 SSS ?! . cosWw, b = 2 SSS f, • cos*pda>, c = 2 SSS f, . cos'ydw, 

h = 2%r l .cos i pcos i ydw i m= 2SS$F,.cos 2 acos 2 7iw, n = 2$$$F,.cos 1 'acos ,! /3rftt>, 

o,= 2SSSF,.cos s acos/3cosy<i<», p i = 2l§v v cos 3 pcosydw, 7 l = 2$$$F,.CGS 3 7COS/3d«, 

a^ =r 2 SSS Fj . cos'o COS7 da>, /3 2 = 2 SSS F r cos ! /3cosa cos7<2«>, y 2 = 2 SSS *V cos 3 7COSadw, 

a 3 :=2SS$sVCOS 3 acos/3e?o), #,= 2$$$F 1 .cos 3 /3cosadw, 7 3 =:2S$SF l .cos !! 7COSacos/fc?». 

It appears, therefore, that the form of the function v for solids (in which we have 
v == v,), is that of a homogeneous function of the second order of the six quan- 
tities, -T-, -tA -; > u > v » w - The most general form of such a function would con- 
tain twenty-one constants, but the function v, which has been just determined, 
contains only fifteen constants, although it has the full number of terms. The 
general equation of equilibrium and motion of solid bodies is therefore 

SSS (*S£ + Yfy + zEf )dm = SSS ly.dxdydz, (12) 

in which the expression (11) for v, must be substituted, and the equation treated 
in a manner analogous to what has been already done for equation (5), this will 
produce the general equations of equilibrium and motion, and also the conditions 
at the limits expressed by means of double integrals. 

It is important to observe that the function v,, which has been just found for 
solids, is quite different from the function v, used by Professor M'Cullagh in his 
Mechanical Theory of Light ; in that theory the function v, which defines the 
medium whose vibrations produce the sensation of light, is a homogeneous func- 
tion of the second order of the three quantities, 
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dy_d£ <%_<% <% _dy 
dz dy dx d£ dy dx* 

this form of the function being deduced from the equation 

d± dv d£ =Q 
dx'dy^dz ' 

which expresses that the medium refuses to change its density. 

It follows from this difference of form of the two functions, that in light the 
molecules of the medium cannot act on each other in the direction of the line 
joining them, for if they did, the function corresponding to the medium should 
contain terms depending on the quantities 

§1 d J. A l v v w 
dx" dy* dz y U ' V ' W ' 

Previously to investigating the differential equations of equilibrium and 
motion, produced by the function v„ it is important to examine whether the 
function v, itself is capable of any simplification in its form by means of a change 
in the direction of the axes of coordinates, to which it is referred ; as, if it be 
possible to refer it to particular axes, for which it becomes simplified, this will 
also simplify the equations of equilibrium and motion depending on it, and will 
prove the existence of three rectangular axes in the body at each point, which will 
have an intimate connexion with the molecular structure of the body at every 
point. The equations for transforming £, ?/, £ are 

£ = a? + brf -f <, 
V = a!% + b'rf + c'£, 

£ = a"r+ b" v '+ &r, 

the six quantities, -jS — t -f, u, v, w, must be transformed by the aid of these 
cuv dy (tz 

three equations, and the transformed values substituted in the function v, ; adopt- 
ing, for brevity, the following expressions, 



p = 2d a", 
p' - 2aa", 
p"=. 2aa', 


q =2b'b", 
q' =2bb", 
q" = 2W, 


r =2c'c" 
r> =2cc", 
r" - 2cc', 


VOL XXI. 


z 





(13) 



160 The Rev. S. Hauohton on the Equilibrium and 

I =b'c"+b"c> t m = c'a" + c"a', n = a'b" + a"b', 
H = b"c + be", rri = c"a + ca", ri = a"b + ab", 
I" = be' + b'c, m" = ca' + c'a, n" = ab' + a'b, 

and, effecting the transformations, we shall obtain the following values for the six 

quantities, -3-, -7-, -~, u, v, w, expressed as linear functions of the six transformed 

d% drf dt , , , 
,».nt,t„s, ^ ^ £,, t, , » , W. 

I = *" » + »'f + <"S + * V '°' + ^ + ° W - 

I =<■"'§? +*"$+ °" 2? + yv '°'+ «'V'.V+ 'W. 

dg' , drf , dt , , , 

' =^<+4£ + 4 + '■>'+ "• v + " v - 

The substitution of these values in the function v, will transform it into another 
form, which still retains the twenty-one terms, though now expressed in terms of 
the six new variables ; the coefficients of which will be known functions of the 
nine angles which determine the position of the new axes of coordinates with 
respect to the old axes. Our object is to discover relations among these co- 
efficients which will enable us, by satisfying them with real values of the nine 
cosines (a, b, c, a', b' t c', a'\ b" y c"), to reduce the function v, to a more simple 
form, and consequently to diminish the number of constants, which, in general, 
is fifteen. 

Let (>„ a„ j„ n„ 3j, j| 4 , tip 3» } 3 ) denote the transformed values of (a fi fi y„ 
a *» &» Ya« a » ft> y 3 ) > the values of these coefficients, when developed, are, 
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N , = {ktfbc + bo'W + ca m b"c") + 3 (lJa'a" + m/W' + nJW) 

+ (a.Cf/i" +p"l') +WI+PH +7*(K + J>'0) ~ (i^ + MiV + n6V) 

+ a,(/a 2 + jt>ftc) + fi^W + pb'c') + 7l (/a"» + pb"d') 

+ atfa* + p'bc) + ft^o" -j- p'bY) + 7a (J'a" 2 + jpW) 

+ a 3 (l"a>+p"bc) + p 3 (l"a'*+p"b'c') + 73 (/"a" 2 + ;>"&"</') 

i, = (a6 3 c + b6'V + cb" 3 </')+3(iJb'b" + ul'bb" + nJ"»') 
+ 30,(0* + ? 6c) + ft(fl' 2 + y*V) + 7l (/6" 2 + qb"c") 
+ «*(f * 2 + tf'&O + 3ft( W 2 + jW) + 7a (Z'6" 2 + q'b"c") 
+ * 3 (/"6 2 + fbc) + ft(/"&' 2 + fb'c') + 3y 3 (?'b"* + y"6"c") 

y l = (ac 3 6 + bc'W + cc" 3 b") + 3(lIc'c" + mI'cc" + n^'cc 7 ) 
+ Sa^/c 8 + rbc) + j8,(fc' 2 + rb'c') + 7l (/c" 2 + r&"c") 
+ ^(tc* + /6c) + 3j8 2 (ZV 2 + r'&V) + 72 (/'c" 2 + I'M') 
+ ^(fV + r^fc) +/3 3 (rc' 8 + r"ftV) + 3y 3 (l"c"* + r"b"c") 

N 2 = (ao'c + bo'V + ca"V') + 3(Lwia'a" + m»»W + nwi'W) 

4- 3a 1 (ma*+pac) + jS^nja' 2 +i>aV) + 7l (W 2 +pa"c") 

+ ^(m'aH j>'ac) + 3ft(m'a /2 +/»'aV) + 72 (m'o" 2 +/>W) (14) 

+ a3(TO , V+^ // ac) + p 3 (»» , V 2 +p"a'c')+ 3 73 (m"o" s +?"o ,, c") 

n 2 = (a6 2 oc + b&'W + cb m a"c") + 3(Lmb'b" + am'bb" + Ktri'bV) 
+ (a,($'m"+?'W)+ p 3 (q"m +qm")+y 3 (qm'+q'm))-(Lac+Ma'c'+xa"c") 

+ a 1 (mb* + qac) + /3.0&' 2 + qa'd) + y l (mb m + qa"c") 
+ a^m'* 2 + (foe) + p i (m'b" i +qfa'c') + y 2 (m'b m + jW') 
+ ^(m'V+^ac) + p 3 (m"b'*+q"a'c')+ y 3 (m"b"* + q"a"c") 

j 8 = (ac 3 o + bc'V + cc"V) + 3(unc'c" + mto'cc" + Nm"cc') 

+ 3a,(«ic 2 + roc) + A(»ic' 2 + ra'c') + 7l (w*c" 2 + ra"c") 

+ a,(wV + /ac) + 3j8 2 (mV*+r'oV) + 7s (mV' 2 + rW) 

+ a 3 (mV+f / 'ac) + A,(w'V i +r"oV) + »y 3 (m"c" i + A"c") 

z 2 
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N 3 = (ao 3 6 + va' 3 b' + ca" 3 b") + 3(LnaV -f vin'aa" + ura'W) 

+ 3a,(na 2 +pab) + ft(«a' 2 + ^a'6') + 7,(raa" 2 + pa"b") 
+ * 2 («'a 2 + /a&) + 3ft(«'a' 2 + j»'a'£') + 7 2 (w'a" 2 +/a"6") 
+ a 3 (Va 2 + />"a&) + /8 3 («'V 2 + p"a!b') + ^{n"a' n +p"a"b") 

3 3 = (Aft'a + b&'V + cb" 3 a") + 3(lwW + m»W + n»"W) 

+ Sa^rai* + qab) + ft(«4' 2 + qa'V) + ^wJ" 2 -f ? a"6") 
+ azin'b* + q'ab) + 3ft(»'6' 2 + ?W) + 7 2 («'6" 2 + q'a"b") 
+ a 3 (»"6 2 + q"ab) + ft(«"6' 2 + ? 'V6') + 3y 3 («"&" 2 + fa"b") 

j 3 = (ac 2 oJ + Bc'Wb' + cc'V'6") + 3(l«cV' + mw'cc" + n»"cc') 
+ («,(/»" +r'W) +ft(r"ra +rn") +y 3 (rn' +r'n))- (Lab + ua'V + aa"b") 

+ a^nc* + rab) + ft(wc' 2 + ra'V) + 7,(»c" 2 + ra"b") 
+ a,(nV + r'aft) + ft(n'c' 2 + rV6') + 7 2 (nV 2 + r'a"6") 
+ a 3 (»V + r"ab) + ft(«"c' 2 -f- r"a'b') + 7 3 («V' 2 + r"a"o") 

In obtaining the foregoing values of the coefficients, use has been made of the 
following relations among the nine angles : 

a = b'c" - b"c', a' = b"c — be", a" = bc' - b'c, 
b = c'a" - c"a\ V - c"a - ca", b" = ea' - c'a, 
c = a'b" - a"b', c' = a"b - ab", c" = ab' — a'b, 

which may be easily seen to be identically true. 

It can be proved from the values of the nine coefficients which have been 
given, that there exists a real system of rectangular axes at each point of the 
body, for which the following equations are true, 

Ni + 3i + Ji = 0, N, + i 2 +^ = 0, f* 3 + 2, + 3 3 = 0, (15) 

for, we obtain from (14), since, p -\- q -f- r =. 0, p' -J- q' + r' zz 0, &c. &c. 

N, + 3, + J, = (A - L)b C + (B - M)b' C ' + (C - *)b"c" 

+ («i + A + y l )l + (« 2 + ft + y*)l' + (« 3 + ft + y 3 )l", 
N* + 2 2 + li = (a — h)ac + (b — m)oV + (c — ti)a"c" 

+ («. +ft + Yi)»» + K +ft +7 2 K + (a 3 + ft +7 3 )ct", 

N 3 + 3 3 + D 3 = (a - l)o& + (b - u)a'V + (c - n)o /, 6" 

+ («. +ft + 7.)» + (^+ft +7>'+ («.'+ A+7.K- 
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Hence, equations (15) are equivalent to 

(A— h)bc + (B - M)b'c' + (C - s)b" C " 

+ («. + A+Y.y + («, + A + 7,y + (»» + & + 7.)' , = 0. 
(a — l)oc -j- (b — u)a'c' + (c — n)a"c" 

+ («i + A+7.)™ + te+A + 7>'+te+A+Y»K' = 0» (16) 
(a - l)o& + (b - u)a'V + (c - N )a"&" 

+ («« + ft + Y.)« + («, + A + Y>' + («, + A + Y>" = 0. 

It is possible to find a unique system of rectangular coordinates for which these 
equations shall be true ; for assume the ellipsoid, whose equation is 

(a — l)^ 2 -|- (b — m)^ 2 + (c — n)z* 

+ 2(a, + ft + 7i )yz + 2(a, + ft + y 2 )xz + 2(a 3 -f ft + 7s )*jf = 1, 

and transform it by the formulae, 

x = ax' -|- by' -\- cz', 
y = oV + jy + cV, 
z =a"x'+b"y'+c"z', 

the coefficients of the rectangles will become the expressions (16), and, conse- 
quently, for axes of coordinates coinciding in direction with the axes of this 
ellipsoid, the equations (15), (16) will be satisfied. 

The function v„ referred to these axes of coordinates, will have the number 
of its constants diminished by three, and consequently the total number will ulti- 
mately be twelve ; it does not appear possible to banish a greater number by any 
transformation with rectangular coordinates, as in this transformation we have 
only three independent variables ; with oblique coordinates, however, we should 
have six, and it is quite possible that there should exist in solid bodies an oblique 
system of axes, which would simplify still further the form of the function v,. 
Such a system of axes, as well as the rectangular system just determined, would 
have an intimate connexion with the molecular structure of the body ; though it 
is not easy to see exactly what that connexion would denote. 

It may be useful to notice here a remarkable relation which exists between 
the function v, and a certain surface of the fourth order, whose properties have a 
great similarity to those of the function ; and the connexion between them is 
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such that a knowledge of any properties of the one must throw a corresponding 
light on the other. This surface is, 

(A* 4 -f- B# 4 + C**) + 6(l#V + M*V + N«y) (17) 

+4yz(3a l x*+{}tf+v l z*) +4a*(a 2 * 2 +3j8 2 y +7^) +4^(*^+j3^ s +37^) = 1 

The relation existing between this surface and the function v x is the following : 
that, if this surface be referred to a new system of rectangular coordinates, the 
new coefficients will be expressed by the same functions of the nine angles as the 
coefficients of the function v„ when similarly changed by a transformation of co- 
ordinates; so that the relations (15) may be established among the coefficients 
of the surface (17)> and vice versa, if we knew any geometrical properties which 
would reduce the surface, there would exist corresponding properties of the 
function v,, which would throw light on the nature of the medium represented 
by v r The relation which exists between the function v, and the surface (17) 
may be proved as follows ; the surface may be put into the following form : 

[>(*»)« + *(ff + c(z 2 )*] + [L(2y*)» + M (2aar)« + n(2^) 2 ] 
+ 2[l(/)(z s ) +m(* s )(«*) +N(*)(y»)] +2[«,(2^)(2*«) + A(2^)(2y*) 

+ 7,(ar*)(2yxr)] 

■(2y*)(o I (aO+ft(y»)+ 7l (*»))- 

+ 2 (2a*)(a,(*») +j3,(y')+ 7 ,(^)) = 1, 

.(2^) («,(*•) +ft(y»)+ 7 ,(««)) J 

which is the same function of the six quantities, a?, y 2 , z\ 2yz, 2xz t 2xy, that v, 
is of -=-, -p-, -r, u, v, w. Also, it may be easily seen that 

x* = aV* + by + <?V 2 + hc.{yz!) + oc(2x'z') + a&(2#y), 
/ = a'V* + 6'y s + </V 2 + JV(tyV) + aV(2*V) + a7>'(2xy ), 

«• = a"V 2 + yy + C " V * + i'V / (2^V)+o , V / (2^V)+o // &' , (2y t y'). (18) 
2yz = px" + qy" + r*' 2 + ?(2^V) + m(2*V) + »(2i;y), 

2ar* = ^V* + ^y i + rV 2 + /'(%V) + m'(2#V) + n'(2x'y' ), 

2*y = />"*" + ^y 2 + r'V 2 + V\2y'sf) + 7»"(2.rV) + n"(2x , y'). 

These equations of transformation, having the same coefficients as the equations 
(13), prove, that the transformation of the surface (17) will produce the same 
coefficients as the transformation of the function v,. 



+ 

+ 

+ 
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I shall now proceed to determine the differential equations of equilibrium and 
motion arising from (12), which is the general equation applicable to all solid 
bodies, and common to them and fluids. 

Substituting for v, its value (11), we obtain, by the calculus of variations, 

y&hrtdxdydz = § (a^ + n -£ + m ■£ + a,u + ^v + a 3 w \ndydz 

+M( a 'S +ft S +7i S +Lo+73V+ftw ^ 
+SK c i +M ^ +L l +7,u+72v+73w )^ rfy 
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-Stf*2+*2+^<*5+»S+*3))w 

Hence, finally, equation (12) becomes 

SSSe(x8? + ySj; + z8£)da%k = A — SSS (p,?I + Q$n + R^dsdydz, ( 19) 

where e denotes the density, A the sum of the nine double integrals, which 
belong exclusively to the limits, and p,, q,, b, the coefficients of 8£, 817, B£, in the 
triple integrals, which give the indefinite equations of equilibrium and motion. 
The equations of equilibrium and motion of solid bodies deduced from (19), are 
the following : 

— eX = P„ — 6Y = Q,, — €Z = b„ (20) 
and 

<P£ , <P*i , d 2 t , /ftlX 

€ ^ = £X + P " € ^ = £Y + Q ' € rf = £Z + B '- < 21 > 

These are the most general equations of equilibrium and motion, and p,, q,, r, are 
expressed in their most general form, without making any supposition as to the 
arrangement of the molecules of the body. 

In order to express the conditions at the limits produced by A, it will be 
convenient to make use of the following transformation of the equation (19). 
Let 

dv. d£ , dn , dV , , , 
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dv t dr , d£ , dn , 

dv. d£ , n drj , dr p 

rfv. d£ , „ dn , dt , „ 

>=?; = +&+*$ + *£ + & + ** + «" 

It appears, from inspection of equation (19). that it may be thrown into the 
form, 

SSS e(x&$ + ySij + tiOdxdydz = SS (nS£ + 3^ + r&Odydz 

+ SS (nfy + ■«£ + 58«)<tofe (22) 

-^((| + | + ^> + (| + | + IMi + g + ?»^ ; 

and, consequently, the equations of equilibrium and motion become 

<£» rfy <fc 



and 



eZ + ^ + l» + ^ =0 
' dz dx dy 

tH> — _i_ ^ a- ^ 4- -^ 

d*« . dz , di , d} _, x 
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If o> represent the element of the bounding surface, and (A, fi, v) the angles 
made by the normal with axes of coordinates, and (H, 9 t ¥ ) the particular forces 
acting at the limits of the solid ; the equation of condition to be satisfied by the 
limits will be 

SS(SS£ -f- QSti -f- ng )u> = A ; 
but we have 

A = SS(**cos\ + 3cosj* -f jocose)^ -f $$ (acos/* + "pcosy -J- ^cosX)whi 

+ SSO C 08 " + DCOSX -f ^COS/ti)tt«£ 

because 

dydz = wcosA, dxdz = tucos/u, dxdy = <o cosv, 
therefore 

S =z NfcosX -{- jcos/t -J- jocose 

$ = 3cos/u + ^cosf -j- jcosA (25) 

f = jcosv + jqcosA -f- "?cos/* 

will be the equations of condition to be satisfied at the limiting surface. If the 
limits be fixed, then we shall have A = 0, without any equation of condition ; 
and if the surface be perfectly free, we shall have 

t*cos\ + 3cos/i + ocosv = 
icos/x -}- bcosv + jcosX = 
J cosy + jq cosA -f- ^> cos/a = 

I shall now proceed to the laws of propagation and transmission of waves in 
an elastic crystalline solid, no external forces acting ; and afterwards consider some 
more particular cases of the theory, which will tend to throw light on the general 
method used in this Paper. 

The differential equations of motion to be integrated are : 

e d?- A d^ + "ay + M ^ + \ ai d^ +a2 dxTz + ^dxTy) 
, {py . a d l r) d^ f <Pt) d 2 ri d*ri\ 

+ a w + A ^ + 7 »a? +2 r5iE + N dlty + &djdz) 

, #£ , . d^ , <?{,./ d% , d* f , d*£ \ 
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d* v _ d*n d l r, #* / d\ d' v d* v \ 

€ W-*df + h a^ + "d? + Y>dWz + ha^ + hd^z) (26) 

e ^ = c d4 + m 3? + L 4+ 2 K<>d±y + *§s + * sfc; 

To these equations of motion we must add the equation of continuity, 
de ,(d$ ,dv,d£\_ 

and we shall have four equations to determine the four quantities (£, i h £, e') as 
functions of (x, y, z, t), e being the original density, and e' the variable density 
at any point. The general integral of these equations is hardly to be expected 
in the present state of the integral calculus; but there exists a particular integral, 
which represents the case of plane waves, and which leads to many important re- 
sults. Let 

£ = cosa.y(0), rj = cosp./(<t>), £ — cosy ./(0), $ — Ix + my + nz - vt, 
where (a, p, 7) denote the direction of vibration of molecules, (/, m, re) the direc- 
tion of normal to wave-plane, and v the velocity of waves. The question to be 
determined is the following : given (/, m, n), to determine, by substituting (£,»/,£) 
in the equations of motion, whether it is possible to find real values of (a, /3, y, »), 
for which such a direction of wave-plane is possible. 

Substituting the values of (£, % £) in the equations (26), we obtain the fol- 
lowing equations of condition : 

e» 2 .cosa = p'eosa + h' cos/3 -f- g'cos?, 

ev\ cos/3 = q'cos/3 -f F'cosy + h'coso, (28) 

etf 2 .cosy = it'cosy + G'cosa -f f'cos/3; 

2 a 2 * 
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where 

p' = a/ 2 + nto 2 -f- m»' 4 + 2a l mn -f 2ajn + 2ajm, 

q' = Bm 2 + lk 2 + n/ 2 + 2/3,«j« + 2/3 2 /ra -f 2ft?*», 

r' = c»* + m/* + l«j ! + 27,/nn -j- 2y Jn + 27./m, (29) 

f' = a/ 2 + ftro* + 7,n 2 + 2l»2/i + 27 3 /w + 2pjm, 

g' = a 2 P + j3 2 m s + 7 2 n 2 -f- 27 3 mn+ 2m/» + 2a/m, 

h' = a.f + /3 3 m 2 + 7 3 « 2 + 2pjmn'+ 2aJ,n + 2nlm. 

Now, it is a known property of surfaces of the second order, that if the equation 
of a surface be 

pV + Q'f + rV + 2p'yz + 2g'xz + 2HVy = 1, (30) 

that the equations (28) will determine the directions and magnitudes of the axes 
of the surface :* the directions being (a, /3, 7), which have three sets of values ; 
and also, if (a, b, c) denote the lengths of the axes, they are determined by the 
equation, 

«_L 1 I 

er ~a 2 ' b» c 2 ' 
ew 2 having three values corresponding to the three sets of values of (a, j3, 7). 

Hence we may deduce a geometrical construction for the direction of vibration 
of molecules corresponding to a given direction of wave-plane. 

Construct the six fixed ellipsoids : 

p = a.t 2 + Ny 2 -f- m* 2 -\- 2a,yz + 2a,jcz + 2a 3 #y = 1, 

q = By" + L2 2 + n^ + 2/3,y* -J- 2&r* + 2/^ry = 1, 

r = cz 2 + m# 2 -f i# 2 + 2 7l #* + 2y#z + 2y#y = 1, (29) 

f = o,j*+ fty«+ 7^ + 2i#* + 2 73 ** + 2ft*y = 1, 

g = a^ 2 + ftjf + 7^ 4- 27^ + 2m** 4" 2a,*y = 1, 

H = a,* 2 4- ^4- 7^ + 2/3^* + 2a,** + 2n*? = 1 ; 

and from their common centre draw the normal to the wave-plane, this will pierce 
the surfaces in six points ; let the corresponding radii vectores be (/»,, p tl , /» /;/ , 
r y , r lt , r jU ), with these construct the ellipsoid, 



x% l £_ j_ * . *(y z jl. XZ y X l\ - 

Yin \ 1 11 

Vid. Leroy, pp. 73, 156. 



r/ "11 Y 111 \ 1 11 //// 
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the directions of the axes of this surface are the three possible directions of mole- 
cular vibration for the given wave-plane ; so that, for a given direction of wave- 
plane, three waves are possible, the corresponding vibrations of the molecules of 
the solid being at right angles to each other. The truth of this construction 
appears from the relations, 

'_! /_ l /_ l /_ l /_ 1 /_ 1 

™ — »» Q — o»K — -j— . F — -5, G — -5-, H — — =- , 
P, P„ P,„ i*. **„ **„ 



III 



which are easily proved. Also, since we may suppose all the parts of the vibra- 
ting body contained in a plane parallel to Ix -J- my -j- nz = 0, to be of the same 
density during the progress of the moving plane which contains them, we shall 
have the three velocities of wave-propagation inversely proportional to the three 
axes of the ellipsoid (30).* 

If wave normals were drawn from the common centre of the ellipsoids (29) in 
every possible direction, and the corresponding ellipsoids (30) constructed for 
each direction, and if on each normal three intercepts were measured, recipro- 
cally proportional to the axes of the ellipsoid, the locus of the extremities of these 
intercepts would be the surface of wave-velocity, or locus of feet of perpendiculars 
from centre on tangent planes of wave-surface ; and the surface formed by radii 
vectores, which are the reciprocals of the wave velocities, will be the surface of 
wave-slowness, or the reciprocal polar of the wave-surface, and a knowledge of its 
properties will serve all the purposes of the wave-surface itself; it is analogous to 
the index-surface in the theory of light. The surfaces of wave- velocity and wave- 
slowness may be determined by the following considerations : 

The cubic equation whose roots are the squares of the reciprocals of the axes 
of the ellipsoid (30) is, 
( P '_,) (q'_ 5 ) ( tt '_,) _ p«( P '_*) - g'Xq'-s) - h'V-s) + 2f'g'h' = 0. 

If in this equation we substitute (x 2 -f- y l -+- z"-) for s, and (->-»-) f° r (A w » »)» 
we shall have the equation of the surface, whose radii vectores measure the wave 

* After having obtained this construction, I learned that M. Cauchy had given a similar one 
for a particular case ; but he has not applied his results to their proper object, and has been mis- 
led in his attempt to account for the phenomena of light by his equations : his method, also, is 
quite different from that used in this Paper. 
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velocities ; and if a radius vector be drawn from the centre, it will pierce this sur- 
face in three points, and the lengths of the three radii vectores measure the three 
velocities of wave-propagation possible for the given direction. This surface, 
however, being only the locus of feet of perpendiculars or tangent planes to 
wave-surface, is of very little use, but it enables us to determine the reciprocal 
polar of the wave-surface ; for, if the radii vectores be made the reciprocals of 
their lengths, the direction remaining the same, we shall have for the locus of 
their extremities the surface of wave-slowness. Its equation is, 

( P y-i) (qV-i) ( B y-i)-FV(py-i)-GV(QV-i)- *YW-i) 

or? +2fVhV = 0, 

(p— 1) (q-1) (r— 1)— f 2 (p— 1)— g 2 (q — 1)-h 2 (r— 1) + 2fgh = 0, (31) 

(p, q, b, f, g, h) being the quantities given in (29). 

This surface is of the sixth degree, and has three sheets, corresponding to the 
three velocities of wave-propagation, and determines not merely the laws of pro- 
pagation of plane waves in a solid, but also enables us to give a construction for 
the direction of waves refracted in passing from one solid into another. 

With a point in the surface of separation (supposed plane) as centre, construct 
the two surfaces (31) for both solids; produce the normal to the incident wave 
to meet its own surface, and from the point in which it pierces it, let fall a per- 
pendicular on the separating plane; this perpendicular will pierce the surface 
(31) of the second solid in three points; the lines drawn from these points to 
the centre are the normals of the three refracted waves, and their lengths are 
inversely proportional to the wave velocities. 

The directions of the internal reflected waves may also be easily found by means 
of this surface ; for we have only to produce the perpendicular backwards to 
meet the surface (31) of the first solid ; and the three lines joining the centre 
with the three points of intersection will be the normals to the three plane waves 
reflected back into the first solid. These and all other constructions for the 
direction of reflected and refracted waves may be easily proved by the properties 
of the wave-surface and the reciprocal properties of the surface of wave-slowness, 
which, for such constructions as these, serves all the purposes of the wave-surface 
itself. The whole theory of wave -surfaces, in light, or in elastic solids, is only a 
development of Huygen's construction for uniaxal crystals. 
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It is important to observe that the directions of the wave and ray, being both 
given, (by the radius vector and perpendicular on tangent plane of the surface of 
wave-slowness), completely determine the direction of molecular vibration in any 
given case. If the direction of the wave alone be given, the problem is indeter- 
minate ; for three parallel tangent planes may be drawn to the wave surface, each 
tangent plane being accompanied by its own direction of molecular vibration, and 
the three directions of vibrations being at right angles to each other : but if the 
direction of the ray be also given, there will be nothing indeterminate, for this 
will show which of the three tangent planes we must select, and, consequently, 
which of the three directions of molecular vibration. 

We now come to the conditions which must be satisfied at the limits ; the 
facilities for obtaining which constitute the great advantage of the method of the 
Mecanique Analytique, as compared with any other method in Mechanics. The 
principal cases of limits which occur are the following: 1st. The limits may be 
completely fixed ; 2nd. The body may have particular forces (E, #, Y) applied 
throughout the whole extent of the limiting surface; 3rd. The limits may be 
perfectly free ; 4th. The limits may be the surfaces separating different kinds of 
bodies, so that the vibration passes from one into the other, changing the laws of 
its propagation at the limits, and requiring particular conditions to be satisfied at 
the separating surfaces. 

1st. If the external surface of the body be completely fixed, so that its 
equation is independent of the time, then there are two cases to be considered : 
the molecules constituting the external surface may be absolutely fixed, so that 
no motion whatever is possible at the limits ; or the external surface may be fixed, 
so that motion perpendicular to it is impossible; but there may be motion along 
the surface itself. In the first case we have no condition at the limits ; for the 
condition to be satisfied at the limits being, 

= A = SS((ncos\ + 3C0S/* -f- acosi/)££ -|_ (acos/* -f- ^cos* + 3cosA)Ej; 

+ Ocosi'-f ocos\ -+- ^>cos^)8£ )w (32) 

is identically true, because 

8£ = 0, fyzzO, ££ = 0. 

But in the second case (which occurs, for example, in hydrodynamics), the 
variations at the limiting surfaces are not zero, but only subjected to the restric- 
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tion that they must take place along a given surface. Let the limiting surface 
he 

f{x, y, z) = 0, 

then the condition to be added to A = 0, is 

and, since (cosA, cos/x, cosv) are proportional to ( -jr-, -j-, -j- J, we shall obtain, 
finally, for the conditions at the limits, 

(33) 

which, by eliminating any one of the variations, will give ultimately the two con- 
ditions to be satisfied. It is evident, a priori, that in this case there must be two 
conditions at the limits, as the equation of the external surface is a function of 
two independent variables. 

2nd. If particular forces (S, $, ¥) act at the external surface, then the equa- 
tion of condition at the limits will be 

A = S(ESf + •&, + «£)«* 
or, 

S = ncos\ + 3 cos/* -j- jocosv 

* = acos/i -|- ^cosv -f- 3cos\ (34) 

Y = ^cosv + ocosA -f- ^COS/i 

which are the same as equations (25). 

3rd. If the limiting surface be perfectly free, then (&£, 817, S£ ) are independent 
of each other, and the general condition A = will give the three conditions, 

N cosA + jcos/t -f- ncosy = 

acos/i + "jcosv + jcosX = (35) 

^COSJ* +QCOS\-t-^COS/i =0 

or, if u =f(x, y, z, t) be the external surface, 
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du , du . du 

dy dz dx 

w du , du , , du 
>Tz+*Tx + ^Ty = ° 

4th. This case, in which we have two media in contact, is the most important 
and difficult of all the cases of conditions at limits. In order to discuss it, it is 
necessary to go back to the general equation of motion and equilibrium (12) ; in 
the case at present considered, the function v, must be conceived to be composed 
of two parts (v,\ v,"), v,' belonging to one body, and v," to the other. Let 
f(x,y,z) = represent the surface passing through the molecules common to both 
media in their state of rest. There will then be three general equations of mo- 
tion, belonging respectively to the molecules, for which 

J( x,y, z)>0, f(x,y,z)<0, J{x t y,z) = 0, 

i. e. according as the molecules belong to the first medium, for which v, = v,' ; 
or to the second, for which v, = v," ; or to both, for which case we have 
v, = v,' -f- v,". These three equations are 

S$ - € " (if 8r + w w ' + %^) dxd y dz = ^< dxd ^ ( 3 7) 

(38) 
= SSS l\(dxdydz + SSS lv"dxdydz ; 

the single and double marks denoting that the quantities belong to the first or 
second body respectively. As the laws of propagation of waves and the limiting 
conditions arising from fixed or free surfaces respectively have been already consi- 
voi. xxi. 2 B 
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dered, we have at present nothing to do with equations (3G) or (37), or with the 
triple integrals in (38), but only with the double integrals in equation (38), pro- 
duced, as in (19), by integrating the right hand member by parts, according to 
the rules of the Calculus of Variations. 

It is easy to see that the condition at the limits will be 

A = A' — A" = 0, 

A', A" arising from the two parts of (38), and having the form of A in equation 
(32). This condition, expressed at full length, is evidently equivalent to the 
equations 

But as the motion of the molecules which constitute the surface f(x,y,z) — 
is the same, whichever medium they are supposed to belong to, we shall have 
evidently ^ = £", */ = »/', g = £" ; and, expressing these in the equations just 
found, and equating to zero the coefficients of the arbitrary variations, we shall 
obtain finally the three following conditions, to be satisfied at the surface which 
separates two contiguous bodies. 

<»'-«">£ + (>'->"> ! + (='-=">* = ° 
( y_a.)* ( y_y.)* + (3f _ „)*=.„ ,39) 

As an example of the formulae, we shall suppose the bodies to be in contact 
along an indefinite plane, which we shall take for the plane (x,y). The surface 
u zzj\x,y, z) = becomes for this case z = 0, and consequently we shall have 
to introduce the conditions, 
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d v __ n du __ du __ 
dx~ ' dy~ ' ' dz ~~ 

The equations of condition to be satisfied at the limits become, therefore, 

»'-»" = o, y_v = o, y-y = o, (40) 

i. e. the three quantities .J, ^, <q, retain their values in passing from one medium 
to another. In order to obtain a more geometrical idea of the nature of these 
equations of condition, in the particular case of plane waves, we shall assume 



_ _d i q _d*r u __<**/ __*£ *_^ 

**-<&' *"<&' *-&'' '~S f ^-rf/ 2 ' 3 ""A»» 

where (p, o, r) are the components of a transversal t,, and (yjg", h) of another 
transversal r /fi so that we have 

p = -j^cosX,, q = Tjcasfi.,, r = ^cos^, 
/ = t^cosX,,, g = t„cos/i„ A = t„cos v 

These quantities (t,, r y/ ) may be called the first and second transversals of the 
real molecular vibration which takes place along (a, /3, 7), and is equal to t; and 
they have a striking analogy to the line that Professor M'Cullagh calls a secon- 
dary transversal in his theory of light ; but there are two such transversals in 
the theory of elastic solids, while there is only one in the theory of light. It 
should be remembered that the equations of condition at the limits (32-40) are 
perfectly general, and independent of the particular integral for plane-waves ; 
but the geometrical constructions, which I proceed to investigate, belong only to 
the particular integral. Previously to determining the meaning of equations (40) 
it will be useful to establish a theorem relating to surfaces of the second order. 
Let the surface be 

as* + ntf + cz* -f- 2nyz + 2e^ -f Txxy = 1 , 
then we shall have 

cos(j>,r') ~ 

pr zz Acosacosa'-f- bcos/Scos^ 4-CCOS7COS7' -f d(cosj3cos7' + cos7COs/3 / ) 
cos^r) -J-E(cos7COsa'-}-cosacos7 / ) -j-F(cosacos/3'-fcos/3cosa / ), (41) 

pr 
where (p,p' t r, r') are the perpendiculars let fall on any two tangent planes and 

2 b 2 
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the corresponding radii vectores, and (a, /3, 7, a, ft, •/) the directions of (r, r'). 

rr«i • f cos(»,r') cos(»', r) . , , „ . ^ ,. 

1 his expression for — — — - = ^ — -, in terms of (a, /3, 7, a , p, 7 ), may 

be proved as follows : the tangent plane to the surface is 

(ax -f- vy -f- E?)^ -j- (By -j" Dz + F^)y -+- (C2 + Ear -J- Dy)z' =c 1 ; 

therefore, if (X, fi, v) denote the direction of the perpendicular, we shall have 

cos\ = /w(ACOsa -f- FCOS/3 -f- ECOS7) 

CO&ft = pr(BCOS/3 -f- DCOS7 -|- fcoso) 

cos* =pr(ccosrf + Ecosa 4- dcos/3) 
therefore 

cosAcosa' -J- cos/icosft -j- COSVCO87' c osVcosa-f-cosM'cosft -j- 'cosy 7 cosy 

pr ~ p'r 1 

= ACOSaCOSa' 4" BC0S/3C0S|3' -f- CCOS7COS7' 

+D(cosj8cos7 / +cos7COS/3')+E(cos7COsa / +cosocos7 / )-|-F(cosacosj3'-|-cos|3cosa ); 

which is equivalent to (41). 

If now we substitute for (£, r), £) their values, 

£ = Tcosa.cos^, t) ■=. tcos/J.cos^, £ = tcos7.cos<^, 

<f> — -r-(/.r -\-iny + hz— ytf) 
in the equations 

^=> = a 4+^+*i+^ u + a ' v+NW 
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we shall obtain, by integrating twice with regard to t, the following expressions, 

p = t / cos\ / = I . k^mip f— t^cosX,, =r I . &,sin^ 

q — Tfio&ft, = I.A^sin^ g — t^cos/*,, = I.A„sin</> (42) 

r =r -^cosy, = I.A^sin^ h =. r /t cosv /t =z l.k in &m<f> 

where we have I = r—r- ; and 

h t = A./cOSa -|- N.OTCOS/3 + M.WCOS7 

-f- ^(ntcosY + ncos/3) + a^ncosa -}- /C0&7) + a,( /cos/3 -+• »racosa) 

h t/ = B.TOCOS/3 -j- 1»»C0S7 -f- N./cOSa 

+ ^(mcosy + wcosjS) -f- j8 2 (wcosa + /COS7) + /3 3 (/cos£ -f- wcoso) 

A„,= c.ncosy -f- M.fcosa -f- L.mcos/3 

+ 7,(mcos7 -)- ncosjS) + 7 2 (raco8a + /C0S7) -f- 7,(/cos/3 + mcosa) 

k t — a,.fcosa -{- /3,.mcos/3 -f- 7,.ncos7 

■+- l(«icos7 -f- «cos/3) -f- 7 3 (ncosa -f- &0S7) -f #,( /cos/3 + wcosa) 

& /y = Oj.&osa + /3 2 .mcos/3 + 7.,.mcos7 

-f- 7 3 (mcoS7 + wcos/3) + M(«cosa + /C0S7) + a,(fcos|8 + wcosa) 

lc m = a 3 ./coso -f- /3 3 .»icos/3 -f- 73.WCOS7 

+ fi^mcosy -f- »cos/3) -f- o,(ncoso -f /C0S7) -|- n(/cos/3 -f rocosa) 

We have, therefore, the following equations to determine the direction and mag- 
nitude of the two transversals : 

k„C08V, — h^COBft, = K C0Sy i/ ~~ */// C0S ^// = ° 

h w cos\ — h t cosv, = ft /w cos\ y/ — Ar^os^ = (43) 

hjcosp, — h„cos\, = k,cosn /t — A^cosX,, = 



T/ = I. Vh] + *J, + h\„ sin^ t„ = I. •*; + AJ, + *?„ sin0 

If we suppose the six fixed ellipsoids (29) constructed, with any point in the 
plane of separation for common centre; and make (p^p^p,^ ^^1,^111) tne P er - 
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pendiculars on tangent planes drawn at the six points where the normal of the 
incident wave-plane pierces the ellipsoids ; the intercepts of the normal itself 
being (r fl r t(t r ul , />,, p lt , p tll ) ; and let (a) denote the direction (a, /3, 7) in which 
the real molecular vibration takes place ; then we shall have the following geo- 
metrical meanings for the quantities (h, &c, k, &c), 



h _ cos(ap,) 
Pf, 


j. _ CM (a*,) 
' *iP, 


k __cos(ap„) 
P„r„ 


k ^COSCOTT,,) 

m ,lpll 


_cos(ap /// ) 
Fur 111 


_cos(a* / „) 



(44) 



Hence the equations for determining the direction of the two transversals become 
finally, 

P,/, r „,- cos («i>„) *»»,— P„ r „- c °< a Pin) C0S M, = 
p,r,. cos(ap,„) cos\,-j> /// r /// .cos(ap / )cosv / = 

p // r ll .cos(ap l )cosfi l — p t r f .cos(ap /l )cos\ = 

(45) 
' r ///P///- cos ( a,r //) cos ''//- ,r /y/ , //- cos ( a ' r //,)cos^ // = 
*, Pr cos («*,„) cos\„ — w,„p„,.coB(air,)cos»'„ = 
*„P„ • cos (aw,) cos?,, — •*,/>, . cos (a-*,,) cos\„ = 

The two transversals are thus completely determined in magnitude and direction, 
and are to be conceived as accompanying the real vibration in its progress, de- 
pending upon its direction and velocity ; and the equations of condition (40) at 
the limits are expressed by saying that, in passing from one medium into another, 
the first transversal preserves its value perpendicular to the separating surface, 
while the second transversal preserves its value parallel to the same surface : 
the real vibration preserving its value both perpendicular and parallel to the 
surface of separation, since the molecules composing it may be viewed as be- 
longing indifferently to either medium. It is important to observe that we 
have thus as many conditions as unknown quantities ; for, the incident vibra- 
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tion being given, we have to determine three refracted and three reflected 
vibrations ; the directions of these are known by means of the surface of wave- 
slowness, and the quantities which are to be determined are the six intensities of 
the reflected and refracted vibrations. Now we have six conditions to determine 
these six quantities ; these conditions are — the conservation of the real vibration 
in passing from one medium into another — the conservation of the perpendicular 
component of the first transversal — and the conservation of the second transversal 
parallel to the plane of separation : the equations of condition given by these 
considerations are, 

r = r, »/ = v", r = r. 

y = y, y = y, »' = »". (46) 

If the vibration pass from a homogeneous solid into one of crystalline struc- 
ture, there will be three refracted and two reflected waves ; and the unknown 
quantities will be the Jive intensities and one direction in the homogeneous body, 
which is not given by the wave-surface ; or, vice versa, if the vibration pass from 
the crystalline into the homogeneous body, there will be three reflected and two 
refracted waves, and the six unknown quantities will be the same as before. Also, 
if two homogeneous bodies be superposed, it will be shown that there are two 
reflected and two refracted waves, and the unknown quantities sxejbur intensities 
and two directions. In fact, whatever be the nature of the superposed bodies, 
the unknown quantities will be six in number, so that they can be completely de- 
termined by the six equations of condition. 

We have hitherto considered the differential equations arising from v, expressed 
in its most general form, and although the relations (15) might have been intro- 
duced, yet the simplification produced by them would be very slight ; the prin- 
cipal effect of these equations being to modify the equation of wave-slowness, and 
it appears very difficult to obtain any other relations among the coefficients of v„ 
which should simplify to any important extent the equations of motion. 

I shall now discuss a few particular cases of solid bodies, in order to show the 
manner in which the general equations already found should be modified when 
applied to such cases. 

I. Let us suppose the solid to have its molecular constitution such, that it is 
symmetrically arranged round three rectangular planes, i. e. the molecular forces 
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are in every respect similar in each of the eight regions into which the three 
planes drawn at any point divide the hody. It is easily seen that in such a me- 
dium the coefficients (a,, /3,, <y„ Oj, $# y 2i a^, fi^ 7 3 ) of the function v, are all zero 
(11), and that consequently v, is reduced to the form, 

^=«D+-Q+«®') + (" ,+ " ,+ " , > 

\ dy' dz dx' dz dx' dy) 

We shall also have the equations (26) reduced to 



e d^ = A a^ + N dJ^ M ^VdsTy + M dxdz) 

#v_# v d> v d> v /<P!; ^^ *e\ 

e d¥-*ay + L ^ + N li* +2 ( L ~ 



( L |fe + ^) < 47 > 



6 



de- dz* +M dx> ±h df + \ M dxdz + L ~dy~dz) 



The equation of the surface of wave-slowness will still be represented by (31), in 
which we must now give to (p, q, b, p, g, h) the following values, 

p = aje* -f- ny* + ms* f = 2iyz 

q == By* + us* -j- N.T 2 g = 2mxz (48) 

r = cz 2 -\- Mar*-f- hy* h == 2N#y 

and its traces on the principal planes of symmetry are given by the following 
equations : 

(i#* + Ma?-1) [(Aar J + Ny— 1) (By* + n**— 1) — 4nV/] = 

(n# 2 + l* 2 — 1 ) [(cz 1 + M^r 2 — 1 ) (a* 3 -f bk»— 1) - 4m**V] = (49) 

(M^ a +n/-i) [( B y -f i^-i) (c^ + l/-i) - 4ityV] = o 

for it appears immediately from the equation of the surface, that its traces are 

(a— 1) [(p— 1) (q-1) — h 2 ] =0 z = Q 
(q-1)[(p-1)(b-1)-g 9 ]=0 y = 
(p— l) [(q-1) (r-1) — f 2 ] =0 x = 
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Hence the traces on the principal planes of the surface of wave-slowness consist 
of an ellipse and a curve of the fourth degree ; which establishes an analogy be- 
tween the case of solid elastic bodies and light, where the corresponding traces 
are a circle and an ellipse. 

The surface of wave-slowness of symmetrical elastic bodies possesses singular 
points, of the nature of nodes, in its principal planes, for the equation of the sur- 
face may be put under any of the forms 

ft +*'•*, = o 
ft + y.+. = o 

ft + *r\+. = 
where the values of 0,, <f> v <f> 3 are 

ft = (H-l)[(P-l)(Q-l)-H*] 
a = (Q-l)[(p-l)(B-l)-G'] 

3 =(p_1)[(q_1)(b-1)-f'] 

Now, if the equation of a surface be 

u = <f> + z\ + = 0, 

it may be shown that, if the surface = has any singular points in the plane 
z = 0, that these are also singular points in the surface u = ; for the conditions 
for singular points are 

du _ ^H — n jg _n. 

dx ~ ' dy ~~ dz~ 

and if z = 0, the part of u of the form (z*.+) will disappear from the three 
equations, 

du du du 

dx dy dz 

which will then be of the form 

dx dy dz 

which are the conditions for singular points in the surface <f> — 0. Hence, if any 
singular points of the surface = exist in the plane z =s 0, these will be also 
vol. xxi. 2 c 
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singular points of the surface u — 0. Applying this principle to the three forms 
of the surface of wave-slowness just given, it will appear that the intersections of 

e— 1=0, (p — 1) (q— 1) — h* = 0, z = 0, 
will be singular points of the surface of wave-slowness in the plane (x, y) ; also 

q— 1=0, (p-1) (b — 1) — g 2 = 0, y = 0, 
will give singular points in the plane {$, z) ; and 

p — 1=0, (q— 1) (b — 1) — p* = 0, x=0, 

will give singular points in the plane (j/, z). 

As the traces in the principal planes consist of a curve of the second and one 
of the fourth degree, we shall have in general eight points of intersection, real or 
imaginary, and therefore the surface of wave-slowness should have twenty-four 
singular points ; but as it is only the real points of intersection that produce any 
effect in the physical problem, it becomes of importance to ascertain the number 
and position of the real singular points. 

I shall first prove that the curve of the fourth degree consists of two ovals, 
lying one inside the other, and not having any point in common. The equation 
of the curve on the plane (x,y) is 

( KX* + vf — 1 ) (htf + Ntf 2 — 1 ) — 4n V/ = 0, 
and its polar equation is 

( Acos 2 a -f- Nsin'a 5 j ( Bsin 2 a -f- ncos 2 o — -» ) — 4N 2 sin 2 acos 2 a = 0, 

which is a quadratic equation with respect to -j; the condition necessary for equal 
roots is 

[(a -f- n) cos 2 a -J- (b + n) sinVj* 
= 4(Acos*a + Nsin 2 o) (usin 2 a -j- ncos 2 «) — 16N 2 sin 2 acqs 2 a ; 

and, if two radii vectores of the curve coincide, i. e. if the curve have a double 
point, this equation of condition must give a real value for a; arranging it with 
respect to tan a, it becomes 
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(b — n) 8 . tan 4 o -f [2<a + n) (b -f- n) — 4ab + 12N a ]tan 2 a + (a - n) 2 = ; 

or, assuming 

u> = (a — n) (b — n) — 4n 2 , 

(b — n) 2 . tan 4 a — 2(o> — 4n 2 ) . taa 2 a + (a — n) 8 = 0. 

Now, in order that this equation should give real values for tana, it must give 
real and positive values for tan*o ; but it can be shown that the roots of this qua- 
dratic equation, if real, are negative, and, consequently, that no real value exists 
for tana. 

The equation solved for tan*a may be put under the form 

tanV(B — n) 2 = (w — 4n 2 ) ± V{* — 4n 2 ) 2 — (a, + 4n*)*, 

or 

tan 2 a. (b - n) 2 = (« - 4n 2 ) ± ^316^. w. 

These equations show that the condition for real values of tan 2 a is, that w must 
be negative or zero ; and that in either case tan 2 a must be negative. Hence we 
have proved, in general, that the two branches of the curve cannot have a real 
point of intersection. The same thing is true of the curves of the fourth degree 
in the other principal planes. 

The geometrical meaning of the condition w = is, that the curve of the 
fourth degree, in each plane, should be the product of an ellipse and circle whose 
equations are 

ax* + By 1 — 1 = 0, *a? -J- Ny 2 — 1 = 0, in plane (x,y). 

B y* -|- as 8 — 1 = 0, Lf + lz 2 — 1 = 0, in plane (y, z). 

cz 2 -f- ax 2 — 1=0, M£ 2 + Ma? 2 —1=0, in plane (jt, z). 

The curve of the fourth degree consists of two branches, each of which cuts 
at right angles, the axes of coordinates, in four points, and the semiaxes of one 

branch are —r=, —-=, and the semiaxes of the other branch are equal, and each 
v a v b 

-/=, and similarly for the other coordinate planes. 

The equation of the ellipse in the plane {x,y) being 

hf + M-r 2 - 1 = 0, 
2c2 
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this curve will evidently cut the oval, whose semiaxes are ~r=, if N be intermediate 

in value between l and m ; supposing, therefore, that m is the mean of the three 
quantities, l, m, n, the ellipse in the plane (x, z) will cut the branch of the oval, 

whose semiaxes are —p= ; and in the other coordinate planes it will lie, in one 
v M 

case, completely outside the oval ; and, in the other case, completely inside the 
oval. Hence there are always at least four singular points on the surface of 
wave-slowness ; whether there be more singular points will depend on the rela- 
tive magnitude of a, b, c, compared with l, m, n; if we assume (as seems probable 
from its being true in homogeneous solids) that a, b, c are greater than l, m, n, 

then the ovals whose semiaxes are ■—=, -7=, -y- will lie completely inside both 

V A * B » C 

the ellipses and the other ovals, whose semiaxes are equal to -7=, — -=, —7=- The 
surface of wave-slowness will therefore consist of three sheets ; one, whose semiaxes 
are -7=, -7=, —7=, isolated, and lying inside the other two sheets ; and the other 

V A V B V C " 

two sheets, having four points in common, like Fresnel's wave-surface, and piercing 

the axis of x in the points -7=, -j=, the axis of y in —7=* —7=, and axis of z in 
r v M v n * Vu v N 

J. 1_ 

v 7 !' vV 

I shall now proceed to discuss some particular cases of the motion, and to 
explain the phenomena which arise from the existence of nodes in the surface of 
wave-slowness. 

Let a plane wave move through the body, perpendicular to z ; the auxiliary 
ellipsoid becomes, in this case, 

m:f* -f- xy 1 + cz ? = 1, 

because /=0, mzzO, n zz 1 ; hence the three possible vibrations of the mole- 
cules will be parallel to the axes of symmetry, and two of them will lie in the 
wave-plane at right angles to each other, and the other will be perpendicular to 
the wave-plane ; and the three velocities of propagation will be proportional to 
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\/c> Vl» -v^m? ^c belonging to the wave whose vibrations are normal to its 
plane, and V\» Vj&, belonging to the waves in which the vibrations are in the 
plane of the wave. Similar conclusions will, of course, hold for waves parallel 
to the other planes of symmetry. It may be observed that, in the case of 
these waves, the direction in which the motion is propagated through the 
solid (which corresponds to the ray in the theory of light) coincides with the 
perpendicular on wave plane ; for it may be proved that the three axes of co- 
ordinates pierce the surface of wave-slowness in eighteen points, at each of which 
the tangent plane is perpendicular to the radius vector ; for the equation of the 
tangent plane is 

« = being the equation of the surface ; but, if we make x = 0, y zz 0, it is 

easily seen that -r- = 0, -^ — 0, and consequently the tangent planes at the 

points where z pierces the surface are given by the equation 

du , , N 
£ (,-,)=0. 

which denotes a plane parallel to (x,y)- 

In general, when a wave-plane passes in any direction through the solid, the 
vibrations of the molecules for which such a direction of wave-plane is possible 
will not lie in the wave-plane ; but there are some cases besides the case of waves 
parallel to the principal planes, in which the vibrations are, one normal to the 
wave-plane, and the other two lying in the wave-plane. 

Let us consider the case of a wave-plane parallel only to one of the principal 
axes (z) ; the auxiliary ellipsoid becomes for this case 

P,* 2 + QuV* + R i^ + 2l W = h 
where 

p, = a/* -\- nm\ ii, = 2slm, 

q, = bwi s -f- til\ because n = 0, 

b, = mP + lw j . 
This equation proves that for every position of waves parallel to an axis, there is 
one vibration in the wave-plane parallel to the axis. 
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Now, the equation which determines the axes of the ellipsoid lying in the 
plane (x, y) is 

H,tan 2 -f- (p, — Q,) . tan0 — H, = 0. 

m 
If we make tan0 =: y, we shall have the condition necessary to be fulfilled, in 

order that the other two vibrations may lie, one in the wave-plane, and the other 
normal to it ; this condition is 

/m[(A — 3n)/ 2 — (b — 3N)m 2 ] = 0. 

The factors / = 0, or m = 0, give the axes ( x, y), which have been already dis- 
cussed, and the third factor gives the equation 

s a— 3x 
tan * = I33V 

which shows that there are two lines in the plane {x, y), making equal angles 
with the axis of x, such that if a wave-plane pass through the solid perpendicular 
to either of them, the three vibrations possible for that direction of wave-plane 
will be two transversal and one normal. 

Thus the existence of nine directions in the solid, for which the vibrations 
are in the wave-plane and normal to it, has been proved ; viz., the three axes of 
symmetry, and two directions in each plane of symmetry. 

To show the effect produced by the existence of nodes in the surface of wave- 
slowness, it will be necessary to consider a wave in its passage from one solid to 
another ; the construction for the refracted wave is as follows : describe the sur- 
faces of wave-slowness (a, b) for both media, having a common centre in the 
plane which separates the two media ; produce the normal to the incident wave 
to meet the corresponding sheet of its own surface (a), and from the point of 
intersection let fall a perpendicular on the separating plane, this perpendicular 
will in general pierce the surface (b) in three points; and the corresponding 
radii vectores will be the normals to the three refracted waves, the perpendicu- 
lars on the corresponding tangent planes being the directions of the refracted 
rays. Suppose the perpendicular should pierce the surface (b) in a node, the 
line joining the centre of (b) with the node will be the normal to refracted 
wave ; but there will be an infinite number of rays, which will form the sides of 
a cone of the second degree, having its vertex at the centre of the surface of 
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wave-slowness, and reciprocal to the tangent cone at the node. Also, there may 
be one ray, and an infinite number of waves ; for if we consider that there are 
four tangent planes of the surface (b), which touch the surface along an ellipse, 
it is evident that there might be a cone of refracted wave normals, of the second 
degree, whose base is one of the ellipses of contact, while there would only be one 
ray, whose direction is the perpendicular on tangent plane : to find, in this case, 
the cone of incident waves which will be refracted into a single ray, we must pro- 
ject, by perpendiculars to the plane of separation, the ellipse of contact of surface 
(b) upon (a) ; then the cone whose base is the projection will be the cone of in- 
cident wave normals ; while the cone whose base is the ellipse will be the cone of 
refracted wave normals, and the unique direction of the ray will be the perpendicu- 
lar to the plane of the ellipse. 

II. If the molecules o( the body be arranged symmetrically round one axis, 
the differential equations of motion will still be the same as (47), with the follow- 
ing relations among the constants : 

a = b = 3n, l = M. 

In these equations the axis ofz is the axis of symmetry; two of these relations 
are evident, a = b, l = m ; and the third may be thus proved. The function 
(f,), in the coefficients (11), is in general a function of (p, 0, <p) ; but since, in 
the present case, everything is similar for all planes passing through the axis of 
symmetry, f, will not be a function of 0, therefore the expressions for a, b, n 
may be integrated once with respect to <p, 

A = B = 2$$$F,COS'Waf, N = 2$$$F,cos 2 acos 2 /3rf« ; 

or, substituting for cosa, cos/3, their values, and assuming 

n = 2§F l sm !i 6p*dpdO, 
we shall have 

A zz b =r \ £l.cos 4 <pd<p, n = \ Xl.sin 2 0cos 2 0<ty>, 

where n is independent of <p ; and, finally, integrating these two expressions, we 

shall obtain 

A = B = 3n. 

These relations introduce corresponding simplifications into the surface of wave- 
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slowness, which will evidently become a surface of revolution round the axis of z. 
There will be no conical refraction for any direction of wave-transmission ; but, 
since two sheets of the surface touch each other on the axis of z, there will be 
only two velocities of waves parallel to (a?, y) t while there will be three velocities 
in every other direction. 

III. If the medium be perfectly homogeneous and uncrystalline we shall still 
use the equations (47) with the additional relations, 

A = B = C, L = M = N, A = 3N, 

which reduce them to the following form : 

df~ \ dx* "*" \df "*" d*P \dxdy "*" dxdzJJ 

d^__ ( *n , (*i d> v \ / d% J£ \\ ( 

<«* v «fa* U*» <¥/ W«fe ^ <&<**// 

These are the differential equations of motion of a solid homogeneous and un- 
crystalline body, and were given first by Navier in the memoir already alluded to. 
Referring to the general form of v„ it will readily appear that, for the case at 
present considered, it becomes 

*.=•(«(£+£+£)+(*•+'■+-> 

4.2^^4. * ^4- * ^Y\ (n\ 

"i" \dy'dz ~*~ dx 'dz T dx' dyjf ^ bV > 

which will give the corresponding equations, 
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which will be found useful in discussing the conditions at the limits for this par- 
ticular case. Previously, however, to examining the limits, it is necessary to 
consider the laws of propagation ; and here I shall confine my attention to the 
particular integral for plane waves already used. In the present case equations 
(29) become 

p = n(3#* +f + z% p' = v(7P -f 1 ), 

q = n(V + *■ + *■). Q' = N(2m 2 +l), 

e = n(3z* + * + f), r' = N(2n* + 1 ), (53) 

f = 2syz, f' = 2n«i«, 

g = 2nxz, g' = 2n/«, 

h = 2tixy, h' = 2n/ot. 

Substituting these values in (28), we obtain 

cw'cosa zr n((2/* + l)cosa -j- 2/mcos/3 -f- 2/«cosy), 

ee 9 cos/3 = N((2m 2 + l)cos/3 -f- 2»mcosY + 2/»»cosa), (54) 

eu 2 cos7 = N ( ( 2» 2 + 1 ) cosy + 2ln cosa -|- 2mn cos/3) ; 

from which we deduce, by multiplying by cosa, cos/3, C0S7, respectively, and 
adding 

e» 2 = n(1 + 2(/cosa + mcos/3 -J- »cosy) 2 ). (55) 

It is evident, on inspection, that cosa = /, cos/3 = m, cosy = w » w '" satisfy equa- 
tions (54), and consequently that one set of vibrations are normal to the wave- 
plane, and the other two in the wave-plane. The equation (55) will give for the 
two cases respectively, 

eo 1 = 3n, 
et) 2 = N, 

the greater velocity belonging to the normal vibration. 

These results may also be deduced from the surface of wave-slowness (31), 
for, by substituting in this equation the values (53), and expanding, it becomes 

[3n(^ +f + z*) - 1] [*(x* +tf + *») - 1]* = 0, (56) 

which shows that it is composed of two concentric spheres, whose radii are 

-7=, -7=, which, therefore, represent the wave-slownesses of the corresponding 

vibrations. 

vol. xxi. 2 D 
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Let the case considered be that of two indefinitely extended homogeneous 
and uncrystalline solids, separated by a plane (which is assumed the plane x, y), 
then the six conditions to be satisfied at the limits are, as in (46), 

r=r'. v' = v", r = r, 

which become, by means of equations (52), 

r=r, v=v", r = r, 

\ dz dx dy ) ~ \ dz dx dy / 

(57) 

The incident vibration may be either transversal or normal ; but as the 
transverse vibration may be resolved into two vibrations, one perpendicular to, 
and the other in the plane of incidence, we must finally, in order to include all 
the possible varieties, discuss the three following cases : 

1. Incident vibration transversal, and perpendicular to plane of incidence. 

2. Incident vibration transversal, and in the plane of incidence. 

3. Incident vibration normal, and in the plane of incidence. 

The plane of separation is the plane of x, y ; z being perpendicular to this 
plane, and the intersection of the plane of incidence with the separating plane 
being the axis of x. 

1st. The incident vibration transversal, and perpendicular to the plane of 
incidence. In this case we shall have £ = 0, t — ; and since the vibration is 
similar along the front of the wave-plane, we shall also have »/ independent ofy. 
The differential equations (50) will therefore be reduced to 

€ ^- N U r+ ^y' w-* VdV+-d^)> (58) 

and the equations (52) will become 
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„ = o, > = ,& 

3=0, »=0, (59) 

, = o, , = .£, 

from which it is evident that the equations of condition at the limits (57) will be 
reduced to the following, 

< = <■• «■% = «■%■• <»> 

recollecting that the condition, # = 0, is implicitly supposed in these latter equa- 
tions, as they belong only to the limits. The equations of condition at the limits 
are thus reduced to two, which is exactly the number required ; for as all the mo- 
tion is perpendicular to the plane of incidence, there will be no normal reflected 
or refracted wave, but only one reflected and one refracted wave, whose vibrations 
are transversal and perpendicular to the plane of incidence, and whose intensities 
are the two unknown quantities in the problem. The equations (59) also prove 
that thejirst transversal is zero, and that the second lies in the plane of inci- 
dence. 

As there are two waves in the first body, and one in the second, we shall have 

r{ = tcos0 + t,cos0,, V' = t'cos^' ; 

t, r y , t being the incident, reflected, and refracted vibration respectively, and 

<j>zz — (lx + nz—vi), <j> t — ^-{1,^ + np—vjk), tf>' = -£ (I'x -f- n'z — v't); 

/, A, v belonging to the incident wave, l fl \,, v, to the reflected, and I', X', v' to 
the refracted. 

Substituting these values in (59), and proceeding as in (42), we obtain 

2ir , ,/m . m . \ 2?r „, „Hri • ,,\ 

~¥'f ~ \\ sm(f> + ir sm v ~¥'f = N \v smtf> ) 

2^ ,, ,/t/ . , rl . , \ 2w .„ „/t'/' . A 

^-. h' = N'^- sin0 + -^ sm 4>,y — . h" = h" [^ sintf>' J ; 

, , ^ , . I X v sinfl , 

but since/ = 4,, vzzv„ XzzX., » = — n., and since r- =—== — == -^-r;, where 

2 D 2 
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0, & are the angles of incidence and refraction, we shall have, when z = 0, 
(ft = y = <f/. Substituting therefore the results obtained above in the equations 
of condition, n' = n", f =/"> we shall deduce finally, 

t + t, = t', 

// \ n n r ' n> 

h'(t-t # ) - = N // .- 7 -j 

or, writing n" = kn', and substituting their values for n, w', \, X', 

T + t, = t', 

T - T ' = T - K ta^ ; < 61 ^ 

whence We obtain 

. ,(. , tanfl\ _ ,/, tanflN 

2t = t(1+k — -), 2t— t'(1 — k — -]; 
\ tanfly ' \ tanfly 

and finally, for the two unknown vibrations, 

__ tanfl' — Ktanfl , _ 2tanfl' _ 

T ' ~ T 'tan6K + Ktanfl ' T ~~ T ' tanA 7 + Ktanfl' *" ' 

To determine the actual position and magnitude of the second transversal, we 

have evidently, for a single wave, 

2 ir „ rasing 

_ /=NT .__ 

2w . hind* 

_A = NT ._; 

whence we obtain 

J = ten0 » ^' ^/M 1 ^ = n .£ sin0, 

i. e. the second transversal lies in the plane x, z, or plane of incidence, and makes 
with the axis of x an angle equal to the angle of incidence. 

2nd. Incident vibration transversal, and in the plane of incidence. In this 
case all the motion will continue in the plane of incidence ; we shall, therefore, 
have rj = ; and reasoning similar to that used in the former case will show 
that | and £ are independent of y. The differential equations of motion will 
therefore become 



Motion of Solid and Fluid Bodies. 1 95 

6 ^- N i 3 ^+^ + 2 rf^> € -df- V \^ + -dF + 2 didz) 

(63) 

1 *- Il i* ,+ * ,+ 55/ "^~ V 3 ^ + ^+ 2 ^/ 

Also, since -=- 4- -jS. = 0, in consequence of the vibration being transversal, 
we shall have the equations 









- A* - N V dz + dx) 



d'h n 



which prove that the Jirst transversal lies in the place of incidence, and that the 
second is perpendicular to it; and also that the equations of condition at the limits 
are the Jour following : 

<4+h=4i*+ d s) <«> 

which are exactly equal in number to the unknown quantities, which are the four 
intensities of the two reflected and two refracted vibrations ; one reflected and 
one refracted wave having its vibration normal to the wave plane, and the other 
two waves having their vibrations transversal. If the incident vibration had been 
normal, it is easy to see that the equations of motion would be the same as 
in the present case, and also that the equations of condition at the limits (65) 
would remain unaltered ; but, in the case of normal vibration, the equation 

-5- 4- -r- = does not hold : the equations (64) will therefore be different, and 
dx dz 

become the following : 
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•=»(■£+$ 



•? = o, 






-KI+I) («) 



3 = 0. 



As the conditions at the limits express that the real vibration is preserved, 
together with the perpendicular component of theirs* transversal and the hori- 
zontal component of the second transversal ; and as two of the four reflected and 
refracted vibrations are normal, and two transversal, it will be necessary, in order 
to express the conditions at the limits, to calculate the two transversals for both 
kinds of waves. 

Substituting, therefore, in equations (64), for | and £, their values for a single 
wave, 



£ = Tcosacos0, 



£ = tcosycos^ 



2tt 



we obtain 



= — (Ix -f- nz — vt), 



2v 2/cosa . 

-3.0:= NT — — 



2v 2»cos'v . 

-^.r = nt — - — - sin*, 



2tt 



/cosy-}-ncosa 



(67) 



and by a similar process from (66), remembering that now the vibrations are 
normal, 



2tt 1 + 2P 
-j-.o = nt. — ^ 

2ir 1 . 

^. ? =NT.-SH10, 



sin^, 



2-n 



2ln 



^ = NT 'X 



(68) 



2tt 14- 2n 2 . 
-^r - nt. — sin*, 
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If now, as in the former case, we designate by (#, <f> f , <p', $ tl , $") the phases of 
the incident, transverse reflected, transverse refracted, normal reflected, and nor- 
mal refracted waves ; and mark in a similar manner the other quantities in the 
problem, we shall obtain for the equations of condition (65), the following, 
(z = 0, making ^ = 4>, = $' = f /t = $", as before) : 

tcoso -|- T^cosa, + T^COSa,, = r'cosa' + T"cosa", 

TCOS7 4" t^osy, + T „ C0S 7// = t'cos^' -f- t"cosy", 

(69) 
2»cos7 2*1,0087, 1 +2»*,_ / ,2n'cosry' „\+2n" 2 \ 

/cos? + ncosa ^a^+n/josa 2Z // «„_ f^ l'cosr/'-^-n'cosa' | /f 2f"n' 

t - -t-t, - + t„— _K^ -, +r -^ 

or substituting for I, n, &c. ... 0, 7, &c, their values in terms of (0, 0,, 0', //S 0"), 
which are all known quantities, we shall have, finally, 

( T -j_ T/ )cos0 + T„sin0„ = t'cos0' + T"sin0", 
(t — T,)sin0 -f t„cos0„ = r'sinO 1 — t"cos0", 
sin20 l+2cos^ / sin20' „ l+2cog' 0"^ 

(- + T /)-X — T "— ^ K (, -v~ T —v 7 — > (70) 

. cos20 8^20,. / / cos20' , „sin20"\ 

(T _ T() __ _ T/# _a = K ^ -^- + T "— j ; 

n" 
k being equal to -7-. 

These four equations determine completely, in terms of known quantities, 
the four unknown vibrations (t,, t, r jt , t"). 

To determine the actual position and magnitude of the two transversals in the 

present case, we must use the equations (67)> which give for the Jirst transversal 

the equations 

r ■nA , sin20 . 

- = - 1, ^y. Vy + r* = Nt . -^=- an*, 

which show that it is situated in the bisection of the angle between z and — x- 
And for the second transversal, 

— — ff = NT.COS20 smd). 
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The position and magnitude of the transversals which have been just deter- 
mined, belong to the transversals of a wave whose real vibration takes place in its 
own plane ; and the position and magnitude of the transversals belonging to a 
normal vibration might be determined in a manner precisely similar from equa- 
tions (68). 

3rd. Incident vibration normal, and in the plane of incidence. The differen- 
tial equations for this case, and the conditions at the limits, are the equations (63) 
and (65). Representing, therefore, by (i>, $,, $', $ tl> <t>") the phases of the inci- 
dent, reflected normal, refracted normal, reflected transverse, and refracted trans- 
verse waves, the conditions at the limits will become the four following, which 
determine, as before, the four unknown intensities : 

(t + T^sinfl -f- T^cosfl,, = T / sin0 / + t"cos0", 
(t — T,)cose + T„sin0„ = t'cos0 / — -r"sin0", 
. l + 2cos 2 sin20„ /,l-f2cosV „ sin20"\ 

The two transversals which I have made use of in this Paper, in investigating 
the conditions at the limits, are merely conceivable mathematical lines, and have 
no real existence in the physical problem : but they are useful to assist in a geo- 
metrical conception of what really happens at the limits. There is one real 
physical vibration, which is to be conceived as accompanied by two mathema- 
tical vibrations, which depend upon the real one — each of the three vibrations 
giving its own conditions at the limits ; so that the two transversals, although 
only conceptions of the mind, are naturally suggested by the conditions of the 
problem, and serve to convey a graphic representation of it. 



